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We study the magnetic excitations of undoped iron oxypnictides using a three-dimensional Heisen- 
berg model with single-ion anisotropy. Analytic forms of the spin wave dispersion, velocities, and 
structure factor are given. Aside from quantitative comparisons which can be made to inelastic neu- 
tron scattering experiments, we also give qualitative criteria which can distinguish various regimes 
of coupling strength. The magnetization reduction due to quantum zero point fluctuations shows 
clear dependence on the c-axis coupling. 
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The discovery of a new class of superconductors with 
transition temperatures exceeding 55K has spurred new 
hope of developinga unified theory of high temperature 
superconductivity^]], 0, S] Like the cuprate superconduc- 
tors, in the iron pnictide compounds superconductivity 
arises from doping a layered antiferromagnet, giving rise 
to tantalizing similarities in the phase diagrams. How- 
ever, there are also striking differences. For example, 
whereas importance is placed on a single d-orbital per 
Cu site in the cuprates, in the iron-based superconduc- 
tors importance is placed on several d-orbitals per Fe site, 
and the total nominal spin per site may be large. Fur- 
thermore, the parent compounds of the iron-based mate- 
rials are semi-metals, rather than Mott insulators as in 
the cuprate case. In addition, recent experiments have 
shown that the electronic couplings in the iron-based su- 
perconductors are more three-dimensional than in the 
cuprate superconductors. [4|, [5|, Ifj 

Although static magnetism tends not to survive in the 
superconducting state of the iron pnictides, magnetic ex- 
citations have been shown to play an important role in 
the superconducting state. In particular, a resonance 
peak has been associated with superconductivity, sug- 
gesting a further connection with cuprate physics. Be- 
cause of the prominent role of magnetism in these ma- 
terials and the connection of magnetic fluctuations to 
the superconducting state, it is important to understand 
the simpler magnetic excitations which are present in the 
parent compound. In order to understand the magnetic 
excitations, we consider an effective Heisenberg model 
with exchange couplings between the net spin associated 
with each site. The effective Heisenberg model may be 
thought to arise from exchange associated either with 
localized magnetic moments, or associated with the net 
moment arising from an SDW associated with itinerant 
electrons. 

At room temperature, most undoped iron-pnictide su- 
perconductors have a tetragonal paramagnetic phase. 
Upon decreasing temperature, the materials show a 
structural transition from tetragonal to orthorhombic. In 
the 122 materials, a three-dimensional long-range anti- 
ferromagnetic order develops simultaneously. This phe- 




FIG. 1: (Color online) Typical magnetic structure of undoped 
iron-based superconductors. 

nomenology constrains effective exchange constants in 
the Heisenberg model, written as [4| 

H = J\a ^ Si • Sj + Jib ^ Si • Sj + J2 ^ Si • Sj 

hj hj i,3 

+ j c 5>.s,-j s ]T(^) 2 (i) 

hj i 

where J\ a and J\b are the nearest neighbor interactions 
along the a- and b-axes, J 2 is the next nearest neighbor 
interaction within the plane, J c is the interaction along c- 
axis, J s is the single ion anisotropy. These couplings are 
illustrated in Fig. [TJ We use linear spin wave theory to 
study the magnetic excitations and sublattice magneti- 
zation reduction due to quantum zero point fluctuations. 

We use Holstein-Primakoff bosons to rewrite the above 
Hamiltonian as [7] 

H = ^ci+^[ Aka k ak+ ^ 5k ^ a -k+ 5 -k tt ka-k)], (2) 
k 1 

where E C \ = {-J\ a + Jib ~ 2J2 - J c - J S )NS 2 is the 
classical ground state energy. The Hamiltonian can be 
diagonalized by using the Bogoliubov transformation 

6k = cosh^kttk — sinh#kaLk? (3) 
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which leads to the result 

H = E C] + E + Y,u(k)blb k 



(4) 



where ou(k) is the spin wave dispersion and Eo is the 
quantum zero-point energy correction. 
The spin wave dispersion o;(k) is given by 



(5) 



where 



A k = 2S(J la - J lh + 2 j 2 + j a + j z + j lh cos ky), (6) 
Bk = 2S( Jio + 2J2 cos cos fc^ + 2J C cos fc^. (7) 



The quantum zero-point energy is then 
£o = ^(-A k + u,(k)). 



(8) 



The presence or absence of gaps at particular points 
in the Brillouin zone may be used to gain qualitative 
information about the state of the system: 



A(7r,0,7r) = 2Sy/J a (2J la + 4J 2 + J s + 2J C ), 



A(7r,7r,7r) = 2Sy/(-2J lh + 4J 2 + J s )(2J la - 2J lb + J s H 
A(0,0,tt) = A(7T,0,7r). 

For example, there can only be a gap at A(7r,0,7r) if 
single- ion anisotropy is present. Measuring a finite gap 
at this point requires that J s be nonzero. In SrFe2As2, it 
has been shown that single- ion anisotropy is present, al- 
though it is a very weak energy scale, J s w 0.015meV.Q 
Given that J s is a small energy scale, we see that mea- 
suring a gap at A(0,7r,7r) would indicate that there 
is anisotropy in the electronic degrees of freedom, i.e. 
J la 7^ Jib- I n this sense, the value of A(0,7r,7r) may 
be taken as a measure of electron nematicity in the 
system. It was furthermore established in Ref. 4 that 
J c ^> J s . When this is the case, the magnitude of the 
gap at A(7r,7r,7r) may be used to indicate proximity to 
the magnetic phase transition at J\ b = 2J2. Because the 
gaps must be real- valued, we see that the system is con- 
strained to have J lb < Min[Ji a + 4f , 2J 2 + 4f ]. Violation 
of this constraint indicates a change in the ground state. 

In the limit of vanishing single ion anisotropy, the gap 
at (tt, 0, 7r) disappears, and the following spin wave ve- 
locities can be defined 



v x = 2S^(J la + 2J 2 )(Jia + 2J 2 + J c ), (10) 
v y = 2S^(2J 2 - J lb )(J la + 2J 2 + J c ), (11) 
v z = 2S v / J c (J la + 2J 2 + J c ). (12) 

Notice that v y becomes imaginary for J\ b > 2J 2l indi- 
cating a change in the classical ground state configura- 
tion. This is consistent with the stability condition noted 




a 



2n 



(9) 



A(0,7r,7r) = 2S\/(2J la - 2J lb + J s ){-2J lb + 4J 2 + J s + 2J C ), 




2J C ), 



(b) 



FIG. 2: (Color online) Spin- wave dispersion band for the an- 
tiferromagnet shown in Fig. [T] (a) Dispersion for J\ a — Jib, 
which corresponds to J\ a — 25, Jib = 25, J2 = 36, J c = 7, and 
Js = 0.05. (b) Dispersion with J\ a ^> Jib, which corresponds 
to Ji a = 40, Jib -5, J 2 20, J c 5, and J s 0.05. 



above which is imposed by requiring that A(0,7r,7r) re- 
main real- valued. 

First principles calculations of the electronic struc- 
ture have made two broad categories of predictions: (1) 



J 2 > Jia - JlbUM, and ( 2 ) 2 ^2 ^ Jia > Jl6.U3,Ul| 

Recent neutron scattering experiments have been used to 
measure the exchange couplings. However, their results 
are quite different. [3,LL2|,[l3[ Here we provide further pre- 
dictions to aid in distinguishing the two cases. Figures 
2(a)| and |2(b)| show the typical spin wave spectrum for the 



two cases. In case (1), there are two small spin wave gaps 
at both (7r,0,7r) and (0,7r, tt). If the system is twinned, 
two spin gaps may be observed. However, in case (2) the 
large interaction anisotropy pushes the spin wave gap at 
(0, 7r, 7r) up to the high energy, which forming a flat zone 
boundary for case (2), and only one low energy spin wave 
gap is expected. 

The neutron scattering cross section is proportional 
to the dynamic structure factor *S(k, cj).[5] In the linear 
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(a) Ji a = 




(b)Ji a > Jib 

FIG. 3: (Color online) Constant-energy cuts (twinned) of the 
dynamic structure factor S(\&,uj) for (a) J\ a — Jib and (b) 
Jia ^> Jib- The x-axis and y-axis correspond to k x and k y 
respectively with the range (— 1.57T, 1.57r). Interaction param- 
eters are same as Fig. [2] 



spin- wave approximation, the transverse parts contribute 
to the structure factor. By symmetry, we have 




FIG. 4: (Color online) S(cj) for the full Brillouin zone. The 
blue line corresponds to case (1) J\ a — Jib, the red line corre- 
sponds to Jia ^> Jib, and the green line is for an intermediate 
case which has Ji a — 40, Jib — 20, J2 — 30, J c = 5, J s = 0.05. 
Interaction parameters in cases (1) and (2) are same as Fig. 
S (cj) is in units of g 2 /J 2 BS e ff- 



energy. In addition, the band top in case (2) becomes 
flat in a very large portion of the Brillouin zone. (See 
Fig- Hp. 

The integrated structure factor S(u) can also be used 
to distinguish the two cases: 



1 1 Ibz 



dk x dk v dk z S aa (k,Lo)S(to - to(k)), 



(14) 

where a = x,y and BZ means integrate over the full 
magnetic Brillouin zone. Numerical results are presented 
in Fig. [H The most dramatic differences are expected 
in the high energy response. There are two broad peaks 
expected for the J\ a = J15 case. There is a sharp peak 
at high energy for the J\ a ^> J15 case which is caused by 
the very large density of states near the spin wave zone 
boundary. As we see from Fig. [2(b)] the spin wave band 
is flat in a large portion of the Brillouin zone. We also 
show a curve of S(co) for the case of J\ a = 2J15 (Ji a is 
not much bigger than Jib). 

The total moment sum rule for a Heisenberg model 
with spin S is defined as [14| 



r(M) = p(k, w ) = s x% 



Ak — Bk 
2cj(k) 



(13) 

where S e ff is the effective spin on an Fe ion, g is the 
g-factor of iron (~ 2), and n(ou) is the Bose occupation 
factor. 

In Fig. [3] we show intensity plots at constant energy for 
the dynamic structure factor 5(k, cc/), assuming a crystal 
with twinned antiferromagnetic domains. In the presence 
of twinning, two concentric spin wave rings are expected 
at low energy if the neutron scattering resolution is high 
enough for case (1). At high energy, the outer ring in- 
creases quicky for case (1) and can form bright spots as 
the rings merge. In case (2), only one ring will be observ- 
able since the energy gap at (0,7r,7r) goes to very high 



„(k)), 



Mo = ^E/^/" d^S aa (k,co) 



= M x + M y + M z 
= g 2 Li%S(S+l). 



(15) 



The transverse part of this corresponds to M a = 
f f f BZ dk x dk y dk z J ^^"^(k, cj), where a = x,y. We 
get M x = My = 0Mg 2 fi 2 B S eff per Fe for case (1) and 
M x = M y = o.580 2 /4S e/ / per Fe for case (2). 

The measured magnetic moment per iron is typically 
less than one Bohr magneton (/i^), which is much smaller 
than the theoretically expected value of ~ 2. 3/i^ per iron 
site predicted by LDA calculations. [1, [TBI, [l6[ In spin- 
wave theory, both the quantum zero point fluctuations 
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the quantum zero point fluctuation. The results show 
that J c can enhance the long range magnetic ordering 
dramatically, which is consistent with the experimental 
findings about the magnetic moments in 1111 and 122 
materials. 

We thank J. P. Hu, P. Dai, J. Zhao, S. Li and E. 
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by Research Corporation and by NSF Grant No. DMR 
08-04748. 



FIG. 5: (Color online ) Jib/Jia dependence of the reduction 
of the sublattice magnetization due to zero point energy of the 
spin waves. Red ( J c = 5 meV) and purple ( J c = meV) lines 
corresponds to J\ a — Jib, blue ( J c = 5 meV) and green (J c = 
0) lines correspond to Ji a ^> Jib- Interaction parameters are 
same as Fig. [2] 



and thermal fluctuations reduce the expected magnetic 
moment per site. Since the energy scale of iron-based su- 
personductors is much larger than the temperature scale, 
we focus on the the quantum zero point fluctuations. The 
sublattice magnetization reduction A^ antnm is defined 
as 



^^quantum 
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dkx dky dkz -^-k 
^tT^tT^tT u;(k) 
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" 2" 
(16) 

In Fig. [5]we present the numerical results for both cases 
(1) and (2). Note that J c has an important effect on 
the magnetization reduction. Empirically, the 122 mate- 
rials are generally more three-dimensional than the 1111 
materials. [13 Thus the effective magnetic moment for 
122 material is expected to be higher than 1111 mate- 
rials. If J c « bmeV, we have Am quanturn « 0.1. For 
S = 1/2, this gives S e ff ~ 0.4, which is close to the 
magnetic moment found by experiments in the 122 ma- 
terials, which have gS e ff « 0.8/i^. [12] 

In summary, we have used the three-dimensional 
Heisenberg model and spin wave theory to study the mag- 
netic excitations, dynamic structure factor, and sublat- 
tice magnetization for the antiferromagnetic spin state 
found in the undoped iron-based superconductors. The 
anisotropy of exchange couplings within the Fe-As plane 
can cause clear differences in the spin wave spectrum, 
expected scattering intensity at constant energy, and in- 
tegrated dynamic structure factor. These can be used to 
determine the degree of anisotropy of magnetic exchange 
interactions within the Fe-As plane. In addition, we have 
calculated the sublattice magnetization reduction from 
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